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Results are obtained for cylindrical leading edges of proposed transatmospheric vehicles by employing a two- 
dimensional viscous shock-layer code for nonequilibrium gas nows. The accuracy and efficiency of the planar 
code is verified through detailed comparisons with other predictions. This study includes results for 6-deg half- 
angle bodies with nose radii ranging from 0.0 1 to 2.0 ft for both cylindrical^ blunted wedges and spherically 
blunted cones (included for comparison). Some results are presented as a ratio of the noncataiytic to the 
corresponding fully catalytic heating value to illustrate the maximum potential for a heating reduction in 
dissociated nonequilibrium flows. Generally, this ratio and the individual heating rates are smaller for cylin- 
drically blunted wedges with small nose radii as compared to the spherically blunted cones (for the same nose 
radius). Therefore, a larger potential exists for heating reduction in cylindrically blunted as compared with the 
spherically blunted surfaces. However, the results presented at higher altitudes (where the slip effects become 
important) show that the spherically blunted nose gives lower stagnation-point heating due to stronger merged 
shock-layer effects as compared with a cylindrically blunted nose. 
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— coefficients appearing in Eq. (5) 

— skin-friction coefficient, 2 T*l(p*.UZr) 
= heat transfer coefficient, 2 q*/(p*UV) 
= freestream specific heat 

= stretching function 

= 1 + «sh fj* 

= r + fjn sh cos 6 

— freestream Mach number 

= molecular weight of mixture 
= coordinate measured normal to the 
body, n*/R% 

= pressure, p*l(pZU¥) 

= heat-transfer rate, q*/{pZUV) 

= nose radius 

= radius of the body surface r*IR% 

= distance measured along the body 
surface, s*/R% 

= temperature, T*IT* vf 

= u*Jic; r , 

= freestream velocity 
= velocity component tangent to body 
surface, uVU* 

= velocity component normal to body 
surface, v*f(JZ 
= axial body coordinate 
= shock angle defined in Fig. 1 
= angle defined in Fig. 1 
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- recombination coefficient for species i 

- stepsize in r? direction 

- stepsize in £ direction 

= transformed r? coordinate, g(f?) 

-- transformed n coordinate, ntn sh 
-- body angle defined in Fig. 1 

- body curvature, 

- coordinate measured along the body 
= density of mixture, 

= shear stress 
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= boundary layer 
= cylinder 
= fully catalytic wall 

- grid-point location along the body 

- noncataiytic wall 

= grid-point location normal to the body 
= perfect gas 
= reference 

- sphere 
= shock 

= viscous shock layer 
= wall 
= freestream 


Superscripts 

* = dimensional quantity 

= quantity divided by corresponding shock 
value 

' = total differentiation with respect to £ 

Introduction 

H IGH-PERFORMANCE transatmospheric vehicles,' 
which fly at hypersonic speeds for both ascent and entry 
conditions, are characterized by small radius leading edges. 
These geometries, while aerodynamically desirable, result in 
large convective heating rates. The leading edges, in particular 
for the wings, must maintain a known surface geometry and 
should not experience a recession of the surface. Conse- 
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quently, the leading edges will have to be actively cooled to 
ensure that these surfaces remain at an acceptable tempera- 
ture. Since these surfaces can constitute many linear feet, the 
required weight of coolant would be substantial. This in- 
creased coolant weight, unfortunately, may translate to a re- 
duction in deliverable payload weight. Therefore, it is im- 
portant to accurately predict the surface heating at the leading 
edges. 

The aerothermal environment about these small leading 
edges is dominated by viscous nonequilibrium effects. The 
investigation of Ref. 2 demonstrated the conservatism in- 
volved in heating calculations based on equilibrium rather 
than a nonequilibrium assumption. The calculations were per- 
formed for a range of nose radii and freestream conditions. 
These calculations, however, were restricted to axisymmetric 
leading edges only. Presently proposed vehicles have leading 
edges ranging in shapes from three-dimensional to two-di- 
mensional (cylindrical). 

The initial purpose of the present investigation was to ob- 
tain nonequilibrium viscous shock-layer results for the two- 
dimensional flows and compare them with the corresponding 
axisymmetric calculations.' 4 While this work does not allow 
a direct assessment of a proposed leading edge, the modified 
code does give the vehicle designer a qualitative assessment 
of surface heating and subsequent cooling requirements for 
the nonaxisymmetric leading edges. 

This study documents a parametric comparison of non- 
equilibrium surface heating levels for sphere-cone (axisym- 
metric) and cylinder-wedge (two-dimensional) shaped 6-deg 
half-angle bodies for different nose radii using the viscous 
shock-layer (VSL) equations. The results are presented pri- 
marily for relatively large Reynolds number conditions, but 
some are also shown (for a sphere and a cylinder) from a 
Navier-Stokes solution procedure 5 with surface slip at low 
Reynolds number conditions. These results provide a com- 
parison between the planar and axisymmetric stagnation-point 
heating under low-density conditions. Emphasis is placed on 
smaller nose radii bodies and lower altitudes. 


Analysis 

Governing Equations 

The conservation equations employed in the present anal- 
ysis are the steady VSL equations for nonequilibrium flow 
over a two-dimensional body at O-deg angle of attack (see 
Fig. I). The normalized form of these equations in the or- 
thogonal, body-oriented, transformed coordinates may be ob- 
tained from Refs. 4, 6, and 7 by setting the variable to 
zero. The major difference between the two-dimensional flow 
and axisymmetric flow- (/ = 1) lies in the form of the global 
continuity equation, which for a two-dimensional flow is 

d d g d 

^ KhPsh“shP«] + ^7 — f^lPshP^ ~ n shVPsh l4 sbPU] = 0 (1) 



are taken from Ref. 9 and are similar to those of Ref. 10. 
Details for the computation of production terms are provided 
in Ref. 4. 

The thermodynamic properties for specific heat and en- 
thalpy, and transport properties for viscosity and thermal 
conductivity are obtained by using curve fits of Ref. 9. The 
semi-empirical method of Wilke 1 ' is used to calculate the gas- 
mixture viscosity. The Prandtl number is set equal to 0.72, 
whereas the Lewis number is assigned a value of 1.4 with an 
assumption of binary diffusion. 

Method of Solution 

The numerical method used for solving the two-dimensional 
VSL. equations is the same as that employed for axisymmetric 
bodies. ' 41 1 However, the shock standoff distance n sU at a 
specified body station m is obtained from the integration of 
Eq. (1) 



“ , [/r,p sh Vsh - < h p sh w sh 


“ Psh l \h(FLUXW)]„ } 



where 


FLUXW = (pv)u (4) 

which is zero in the absence of injection through the surface. 

Equation (3) is different than that for an axisymmetric flow, 
where a quadratic equation for w sh is obtained. 4 

The first-order global continuity equation [Eq. (1)] is re- 
written (at a given body station m) as 

(<*X , 1,2*V„ , 1 + 4 „ 2 V f „ „ + (tX . I’lPm.rt ♦ « 


where 


+ (4 )„ . 1/2 P mM = (On , 1.2 


h i = l+ n sh rjK 


( 2 ) 


with the help of the equation of state 


( 5 ) 


Boundary Conditions and Flowfield Properties 

Most of the results presented here are with the no slip 
boundary condition at shock and body. The shock boundary 
condition for no slip case is obtained from the Rankine-Hu- 
goniot relations. The shock and surface slip boundary con- 
ditions (implemented for low density flows) are obtained from 
Refs. 7 and 8, respectively. 

For nonequilibrium airflow, a seven-species (N-,, CK, N, O, 
NO, NO 4 , and e ) chemical model is used to evaluate the 
species production terms that appear in the energy and species 
continuity equations. The expressions and the reaction rates 
used to compute the production terms in the present study 


P = (6) 

which is employed to eliminate the density p from Eq. (1). 
Equation (5) is similar in form to that given in Ref. 4 for an 
axisymmetric flow and is solved in a similar coupled way with 
the normal momentum equation. The coefficients appearing 
in Eq. (5), however, are different from those given in Ref. 4 
and are provided in Appendix A of Ref. 12. The coefficients 
for the /7-momentum equation are the same as those given in 
Ref. 4. 

A procedure similar to that of Ref. 4 is also used to obtain 
the stagnation-point solution. This solution is obtained by 
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expanding the various flowfield quantities in terms of distance 
£ along the body surface. These expansions reduce the partial 
differential equations to ordinary differential equations at the 
stagnation-point singularity (£ = 0) which may be integrated 
normal to the surface with appropriate shock and surface 
boundary conditions. The solution procedure 4 then marches 
along the body using the stagnation-point solution. 

The Vigneron condition' is used for the pressure gradient 
in the streamwise momentum equation in the subsonic nose 
region. The global iterations are employed in this region to 
update the pressure (and shock shape). In the supersonic flow 
region, the pressure gradient term is approximated with a 
two-point backward difference, and a predictor-corrector 
approach' is used for the shock shape to march the solution 
downstream. It should be mentioned that unlike the axisym- 
metric calculations, 3 h converged results cannot be obtained 
for the two-dimensional shape without the use of Vigneron 
condition. 


Results and Discussion 

Detailed results from nonequilibrium calculations are pre- 
sented in this section. First, a comparison with existing nu- 
merical predictions is provided. Results for a cylindrical ly 
blunted slender wedge and a spherically blunted cone for a 
range of nose radii are presented next. Some results are also 
provided as a ratio of the noncatalytic heating to the corre- 
sponding fully catalytic value. This ratio (which is also a mea- 
sure of the degree of flowfield nonequilibrium) demonstrates 
the maximum potential for a surface heating-rate reduction 
in the presence of dissociated nonequilibrium flow. Finally, 
some results from the planar and axisymmetric flow calcu- 
lations are also included at higher altitudes using the stag- 
nation region Navier-Stokes code of Ref. 5 with surface slip. 
Reference 12 contains comparisons and other detailed cal- 
culations for the perfect-gas flows. 




rig. X VUIII(lal IMJII Ul |Ti cuiviiu iv.iuiu 

a 6-deg cylindrically blunted wedge: a) surface pressure distribution 
and b) surface heat transfer rate distribution. 

of the predicted distributions due to overexpansion and re- 
compression of the flow over a spherically blunted cone. 


Comparison with Other Predictions 

Present predictions are compared with those obtained from 
LAURA 13 14 in Figs. 2a and 2b for a fully catalytic wail. Both 
the calculations give similar qualitative results. The maximum 
difference between the surface pressure predictions (Fig. 2a) 
is about 10%, whereas for the surface heat transfer rate (Fig. 
2b) it is about 20%. These differences are believed to be due 
to the different transport and thermodynamic properties em- 
ployed in the two predictions. Presently, it is being planned 
to implement the properties from Ref. 9 (employed in the 
present calculations) in code LAURA. 

Comparison with Axisymmetric Calculations 
and Parametric Studies 

Two-Dimensional and Axisymmetric Calculations 

Results from a comparative study of the cylindrically blunted 
wedge and spherically blunted cone are shown in Fig. 3. These 
calculations are for 6-deg half-angle noncatalytic bodies with 
a nose radius of 0.125 ft along a surface distance of 300 nose 
radii. The shock inflection shown in Fig. 3a for the sphere- 
cone calculations is not noticeable for the cylinder-wedge pre- 
dictions. The effect of this shock inflection is seen in the 
surface pressure distribution of Fig. 3b and heating rate dis- 
tribution of Fig. 3c. For cylinder-wedge calculations, the sur- 
face quantities continuously decrease beyond the tangency 
point. In the stagnation region and its immediate vicinity, the 
surface pressures are almost same for the two calculations as 
shown in Fig. 3b. However, due to the larger shock standoff 
distance for the cylindrically blunted wedge as compared to 
a sphere cone (Fig. 3a), the surface heating-rate values are 
generally lower than those predicted for a spherically blunted 
cone. There is a slight crossover (not noticeable in Fig. 3c) 


Bluntness Effects 

Similar to the studies of Refs. 2 and 15, nose bluntness 
effects are studied through the ratio of noncatalytic to fully 
catalytic heating rates. This ratio shows the maximum poten- 
tial for a heating reduction in dissociated nonequilibrium flow , 
with a smaller value of the ratio implying a larger potential 
for reduction. For a finite catalytic wall this ratio will be larger 
and, therefore, the actual heating reduction would be less 
than the values shown in Fig. 4. This figure gives the heating 
reduction ratio through a wetted distance of 30 ft for nose 
radii values of 0.125 and 0.25 ft. The heating-reduction ratio 
generally decreases to a value of less than 0.40 over the spher- 
ical (or cylindrical) nose portions of the body. This ratio stays 
less than this value for a cylinder-wedge for running lengths 
of less than 30 ft. The ratio becomes as large as 0.9 for a 
sphere-cone body. 

It may be noticed from Fig. 4 that the heating ratio for a 
sphere cone rapidly approaches a value of about 0.9 on the 
downstream side for the smaller (R% = 0.125 ft) nose radius. 
For such large values of the ratio, a finite-catalytic wall bound- 
ary condition would have a very small influence on the heat- 
ing. A large value of the heating ratio for a low angle sphere- 
cone is due to the fact that the effects of dissociated species, 
which are produced primarily in the stagnation region, di- 
minish on the surface heating as the flow spreads over a larger 
surface area at these conditions. 

For the two nose radii considered in Fig. 4, the heating- 
rate distributions for a fully catalytic wall are given in Fig. 5. 
This figure when used with Fig. 4 can give the individual 
heating reduction for a given nose radius. The heating rate 
values for a cylinder-wedge, in general, are lower than the 
sphere-cone values (by as much as 50% for R*„ = 0.125 ft 
and s* = 30 ft). For the larger nose radius (R% = 0.25 ft). 
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Fig. 3 Comparison of shock standoff distance and surface quantities 
for 6-deg two-dimensional (cylinder-wedge) and axisymmetric (sphere- 
cone) bodies in nonequilibrium flow: a) shock standoff distance, b) 
pressure distribution, and c) heat transfer rate distribution. 



Sf , ft 


Fig. 4 Nonequilibrium heating ratio vs wetted distance. 



b) S*,ft 

Fig. 5 Nonequilibrium heat-transfer rate distribution for different 
nose radii: a) R% = 0.125 ft and b) = 0.25 ft. 

differences between the two predictions decrease. Reference 
12 also provides results for R% = 0.5 and 1.0 ft. 

Stagnation- Point Heating vs Nose Radius 

The heating-reduction ratio as a function of nose radius is 
shown in Fig. 6a. The heating ratios for both cylinder and 
sphere calculations first decrease to minimum values (maxi- 
mum nonequilibrium conditions), and then increase with the 
increase in nose radius. This trend with increasing radius is 
due to the flowfield chemistry changing from frozen state to 
equilibrium condition. A lower value of the heating ratio 
implies a higher degree of nonequilibrium. For flows with a 
low heating-ratio value, considerable reduction in fully cat- 
alytic wall heating value can be obtained by making the sur- 
face noncatalytic. Figure 6a shows that the flowfield chemistry 
“freezes” slower with the decrease in nose radius for a cyl- 
inder (due to large relaxation distance) as compared with a 
sphere. Therefore, a bigger potential exists for the surface 
heating reduction for smaller nose radii cylindrically blunted 
bodies. 

The heating-reduction ratio of Fig. 6a shows a potential for 
reduction. It does not, however, give actual values for the 
heating reduction. These values are given in Fig. 6b. Figure 
6a gives the same two-dimensional value of about 0.5 for 
heating reduction for nose radii of about 0.01 and 0.06 ft. 
However, the actual noncatalytic wall (NCW) heating value 
in Fig. 6b at R% - 0.01 ft is about 3 times the value at R„ 
= 0.06 ft. Similarly, the same axisymmetric heating-reduction 
ratio of 0.5 for = 0.04 and 0.15 ft (from Fig. 6a) gives 
actual values which differ by a factor of 1.5 (Fig. 6b). Thus, 
Fig. 6b gives the required heating values to a designer. How- 
ever, the fully catalytic wall values of Fig. 6b can be reduced 
by making the surface noncatalytic as shown in this figure. 
The percentage of reduction for a given nose radii can be 
obtained from Fig. 6a. 

Figure 6c compares the stagnation point heating on a spher- 
ical nose for a range of nose radii. The present VSL heating 
rates obtained for equilibrium flow chemistry are almost sim- 
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Fig. 6 Stagnation-point quantities vs nose radius: a) heating ratio, 
b) heating rates, c) heating rates (for sphere), and d) heating ratio. 


ilar to those predicted for nonequilibrium flow chemistry with 
the fully catalytic wall boundary condition, except for the 
smaller (R% < 0.2 ft) nose radii. Clearly, the present heating 
results do not scale as 1 /\ZR% for smaller values of the nose 
radii due to large viscous and nonequilibrium effects. The 
classical boundary-layer result of Fay and Riddell is also 
shown along with the coupled boundary layer and inviscid 
calculations 17 for equilibrium flow chemistry. The coupled 
calculations give results which are within 4 and 17% of the 
equilibrium VSL results at nose radii of 2.0 and 0.02 ft, 
respectively. These calculations are obviously in better agree- 


C H 



3000 2000 1000 600 300 200 100 60 30 20 10 
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Fig. 7 Comparison of stagnation-point heat transfer coefficient for 
a cylinder and a sphere with slip. 


ment with the VSL predictions than the classical boundary- 
layer results of Fay and Riddell. Since the results for the two- 
dimensional (cylindrical stagnation point) calculations scale 
as 1/V2 with the corresponding axisymmetric (spherical nose) 
computations for equilibrium flow chemistry, 18 differences in 
the various heating rates for the two-dimensional computa- 
tions are expected to be similar to those shown in Fig. 6c. 

Present calculations for nonequilibrium flow chemistry with 
noncatalytic and fully catalytic walls are intended to provide 
lower and upper bounds, respectively, for the heating rates 
to a designer. Therefore, no results are presented as a function 
of the recombination rate parameter as provided in Ref. 16. 
In addition, since a corresponding approximate nonequilib- 
rium analysis (for use with the Fay and Riddell’s correlation 
for frozen flow) is not available, a comparison of Fay and 
Riddell’s results with the VSL nonequilibrium calculations is 
not provided. 

The stagnation-point heating for a cylinder is, generally, 
lower than that for a sphere as shown in Fig. 6d. For a non- 
catalytic surface, the heating rate for a cylinder q l is 35% of 
the sphere value q s at R% = 0.01 ft, and increases to about 
75% at R% = 0.5 ft before decreasing slightly. For a fully 
catalytic wall, the cylinder value increases from about 50 to 
70% of the sphere value over the range of nose radius shown 
in Fig. 6d. The cylinder values range between 60-70% of the 
sphere values for the perfect gas calculations 12 shown here. 
Also given in this figure is a value of l/\/2 for the equilibrium 
and perfect-gas flows. 18 For larger nose radius, nonequilib- 
rium values of the qjq s ratio appear to approach this value. 


Stagnation-Point Results at Higher Altitudes 

The viscous shock-layer results presented so far are for 
lower altitude flight conditions and have been obtained with- 
out the slip conditions. Some stagnation-point surface heat 
transfer results have also been obtained at higher altitudes, 
where the slip conditions become important. The present vis- 
cous shock-layer method with surface 8 and shock slip 7 is em- 
ployed at an altitude of 270,000 ft, and a Navier-Stokes so- 
lution procedure 5 with surface slip 8 is used for altitudes equal 
to and greater than 270,0(X) ft. The two methods give almost 
similar values for the heat transfer coefficients for both a 
cylinder and a sphere at the common altitude of 270,000 ft as 
shown in Fig. 7. One interesting point to note in this figure 
is that stagnation-point heat transfer coefficient for a sphere 
becomes smaller than the corresponding value for a cylinder 
for altitudes greater than about 315,000 ft. This is a conse- 
quence of the stronger merged shock-layer effect on a sphere 
as compared with a cylinder. 


Concluding Remarks 

Results are presented from the study of two-dimensional 
nonequilibrium VSL flows. These are among the earliest VSL 
results for such flows reported in the literature. Comparisons 
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of the present results are made with those obtained from the 
Navier-Stokes LAURA algorithm, and the results compare 
well. These comparisons indicate that the present planar code 
is as accurate and requires as few computer resources as its 
axisymmetric counterpart when compared to LAURA. The 
axisymmetric flow results are included in the present study to 
highlight the differences between the planar and axisymmetric 
flows. 

The investigation includes a detailed study of the nose 
bluntness effects on the stagnation point as well as down- 
stream nonequilibrium flow conditions. The downstream in- 
fluence of nonequilibrium flow for a small nose radius two- 
dimensional body is found to be much smaller than for the 
larger nose radii. Although a ratio of noncatalytic to fully 
catalytic heating rates was used to characterize the maximum 
potential for a nonequilibrium heating reduction over the range 
of conditions, a comparison of individual stagnation heating 
rates demonstrated that a comparison of the same ratio value 
was not indicative of the actual heating reduction. Further- 
more, the fully catalytic wall heating rates (which are closer 
to the equilibrium value) illustrate the amount of conservatism 
contained in such a calculation, as compared with a noncat- 
alytic (or finite-catalytic) wall value. This conservatism would 
result in excessive coolant requirement at the cost of deliv- 
erable payload. 

Generally, the heating-rate values for a cylindrically blunted 
wedge are lower than those for a spherically blunted cone, 
except at very high altitudes. Under those low-density con- 
ditions the merged shock-layer effects seem to be stronger for 
a sphere than for a cylinder. This results in lower heating for 
a spherical than for a cylindrical stagnation point. 
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